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Abstract. We propose to use alternative cosmic tracers to measure the dark energy equation 
of state and the matter content of the Universe [w(z) & fl m ]. Our proposed method consists 
of two components: (a) tracing the Hubble relation using Hll-like starburst galaxies, as an 
alternative to supernovae type la, which can be detected up to very large redshifts, z ~ 4, 
and (b) measuring the clustering pattern of X-ray selected AGN at a median redshift of ~ 1. 
Each component of the method can in itself provide interesting constraints on the cosmological 
parameters, especially under our anticipation that we will reduce the corresponding random 
and systematic errors significantly. However, by joining their likelihood functions we will be 
able to put stringent cosmological constraints and break the known degeneracies between the 
dark energy equation of state (whether it is constant or variable) and the matter content of 
the universe and provide a powerful and alternative rute to measure the contribution to the 
global dynamics and the equation of state of dark energy. A preliminary joint analysis of X-ray 
selected AGN (based on a small XMM survey) and the currently largest SNIa sample (Kowalski 
et al 2008), provides: f2 m = 0.28±g;{^ and w = -1.0 ± 0.1. 



1. Introduction 

We live in a very exciting period for our understanding of the Cosmos. Over the past decade the 
accumulation and detailed analyses of high quality cosmological data (eg., supernovae type la, 
CMB temperature fluctuations, galaxy clustering, high-z clusters of galaxies, etc.) have strongly 
suggested that we live in a flat and accelerating universe, which contains at least some sort of 
cold dark matter to explain the clustering of extragalactic sources, and an extra component 
which acts as having a negative pressure, as for example the energy of the vacuum (or in a more 
general setting the so called dark energy), to explain the observed accelerated cosmic expansion 
(eg. Riess, et al. 1998; 2004; 2007, Perlmutter et al. 1999; Spergel et al. 2003, 2007, Tonry 
et al. 2003; Schuecker et al. 2003; Tegmark et al. 2004; Seljak et al. 2004; Allen et al. 2004; 
Basilakos & Plionis 2005; 2006; Blake et al. 2007; Wood-Vasey et al. 2007, Davies et al. 2007; 
Kowalski et al. 2008, etc). 



Due to the absence of a well-motivated fundamental theory, there have been many theoretical 
speculations regarding the nature of the exotic dark energy, on whether it is a cosmological 
constant, a scalar or vector fields which provide a time varying dark-energy equation of state, 
usually parametrized by: 

PQ = w{z)p Q , (1) 
with pq and pq the pressure and density of the exotic dark energy fluid and 

w(z) = w + wif(z) , (2) 

with wq = w(0) and f(z) an increasing function of redshift [ eg., f(z) = z/(l + z)] (see Peebles 
&: Ratra 2003 and references therein, Chevalier &; Polarski 2001, Linder 2003, Dicus & Repko 
2004; Wang & Mukherjee 2006). Of course, the equation of state could be such that w does not 
evolve cosmologically. Current measurements do not allow us to put strong constraints on w, 
with present limits — 0.8 (eg. Tonry et al. 2003; Riess et al. 2004; Sanchez et al. 2006; 
Spergel et al. 2006; Wang & Mukherjee 2006; Davies et al. 2007). 

Two very extensive recent reports have identified dark energy as a top priority for future 
research: "Report of the Dark Energy Task Force (advising DOE, NASA and NSF) by Albrecht 
et al. (2006), and "Report of the ESA/ESO Working Group on Fundamental Cosmology", by 
Peacock et al. (2006). It is clear that one of the most important questions in Cosmology and 
cosmic structure formation is related to the nature of dark energy (as well as whether it is the 
sole interpretation of the observed accelerated expansion of the Universe) and its interpretation 
within a fundamental physical theory. To this end a large number of very expensive experiments 
are planned and are at various stages of development. 

Therefore, the paramount importance of the detection and quantification of dark energy for 
our understanding of the cosmos and for fundamental theories implies that the results of the 
different experiments should not only be scrutinized, but alternative, even higher-risk, methods 
to measure dark energy should be developed and applied as well. 

1.1. Methods to estimate the dark energy equation of state. 

A large variety of different approaches to determine the cosmological parameters exist. A few 
of the most important ones are listed below: 

• CMB power-spectrum + Hubble relation: By measuring the curvature of the universe (Ofc) 
from the CMB angular power spectrum and using SN la as standard candles to trace the 
Hubble relation (eg. Riess et al. 2004) to measure a combination of VL m and then the 
values of all the contributing components to the global dynamics can be constrained (using 
the fact that fi^ + O m + 17a = 1) 

• Hubble relation + Clustering of extragalactic sources: If, however, the dark energy 
contribution is not due to a cosmological constant but rather it evolves cosmologically, 
then the Hubble relation provides a degenerate solution between the contribution to the 
global dynamics of the total mass and of the dark energy, even in the case of an almost flat 
geometry. In order to break this degeneracy one needs to introduce some other cosmological 
test (eg., the clustering properties of galaxies, clusters or AGN, which can be compared with 
the theoretical predictions of an a priori selected power-spectrum of density perturbations 
- say the CDM - to constrain the cosmological parameters such as Q m , w, a% and h; eg., 
Matsubara 2004). 

• Baryonic Acoustic Oscillation method, which was identified by the U.S. Dark Energy Task 
Force as one of the four most promising techniques to measure the properties of the dark 
energy and the one less likely to be limited by systematic uncertainties. BAOs are produced 
by pressure (acoustic) waves in the photon-baryon plasma in the early universe, generated 



by dark matter (DM) overdensities. At the recombination era (z ~ 1100), photons decouple 
from baryons and free stream while the pressure wave stalls. Its frozen scale, which 
constitutes a standard ruler, is equal to the sound horizon length, r s ~ 100 h~ l Mpc 
(e.g. Eisenstein, Hu &: Tegmark 1998). This appears as a small, ~ 10% excess in the 
galaxy, cluster or AGN power spectrum (and 2-point correlation function) at the scale 
corresponding to r s . First evidences of this excess were recently reported in the clustering 
of luminous SDSS red-galaxies (Eisenstein et al. 2005, Padmanabhan et al. 2007). A 
large number of photometric surveys are planned in order to measure dark energy (eg., 
the ESO/VST KIDS project, DES: http://www.darkenergysurvey.org, Pan-STARRS: 
http : //pan-starrs . if a.hawaii . edu). 

• Other important cosmological tests that have been and will be used are based on galaxy 
clusters. For example, (a) the local cluster mass function and its evolution, n(M,z), which 
depends on f2 m , on the linear growth rate of density perturbations and on the normalization 
of the power-spectrum, erg (eg., Schiiecker et al. 2003; Vikhlinin et al. 2003; Newman et 
al. 2002; Rosati et al., 2002 and references therein), (b) the cluster mass-to-light ratio, 
M/L, which can be used to estimate Q m once the mean luminosity density of the Universe 
is known, assuming that mass traces light similarly both inside and outside clusters (see 
Andernach et al. 2005 for a recent application), (c) the baryon fraction in nearby clusters 
(eg., Fabian, 1991; White et al., 1993). Assuming that it does not evolve, as gasdynamical 
simulations indicate (eg. Gottldber & Yepes 2007), then its determination in distant clusters 
can provide a geometrical constraint on dark energy (Allen et al., 2004). 

1.2. Objectives of our Approach 

We wish to constrain the dark energy equation of state using the combination of the Hubble 
relation and Clustering methods, but utilizing alternative cosmic tracers for both of these 
components. 

From one side we wish to trace the Hubble function using HQ- like starburst galaxies, which 
can be observed at higher redshifts than those sampled by current SNIa surveys and thus at 
distances where the Hubble function is more sensitive to the cosmological parameters. The 
HII galaxies can be used as standard candles (Melnick, Terlevich & Terlevich 2000, Melnick 
2003; Siegel et al. 2005) due to the correlation between their velocity dispersion, metallicity 
and Hp luminosity (Melnick 1978, Terlevich & Melnick 1981, Melnick, Terlevich & Moles 1988). 
Furthermore, the use of such an alternative high-z tracer will enable us to check the SNIa 
based results and lift any doubts that arise from the fact that they are the only tracers of the 
Hubble relation used to-date (for possible usage of GRBs see Ghirlanda et al. 2006; Basilakos 
& Perivolaropoulos 2008)0. We therefore plan to improve the L(Hp) — a distance estimator 
by investigating all the parameters that can systematically affect it, like stellar age, metal and 
dust content, environment etc, in order to also determine with a greater accuracy the zero-point 
of the relevant distance indicator. The possibility to use effectively HH-like high-z starburst 
galaxies as cosmological standard candles, relies on our ability to suppress significantly the 
present distance modulus uncertainties {a^ ~ 0.52 mag; Melnick, Terlevich &: Terlevich 2000), 

GRBs appear to be anything but standard candles, having a very wide range of isotropic equivalent luminosities 
and energy outputs. Nevertheless, correlations between various properties of the prompt emission and in some 
cases also the afterglow emission have been used to determine their distances. A serious problem that hampers 
a straight forward use of GRBs as Cosmological probes is the intrinsic faintness of the nearby events, a fact 
which introduces a bias towards low (or high) values of GRB observables and therefore the extrapolation of their 
correlations to low- z events is faced with serious problems. One might also expect a significant evolution of the 
intrinsic properties of GRBs with redshift (also between intermediate and high redshifts) which can be hard to 
disentangle from cosmological effects. Finally, even if a reliable scaling relation can be identified and used, the 
scatter in the resulting luminosity and thus distance modulus is still fairly large. 



which are unacceptable large for precision cosmology. 

Prom the other side we wish to determine the clustering pattern of X-ray selected AGN at a 
median redshift of ~ 1, which is roughly the peak of their redshift distribution (see Basilakos et 
al. 2004; 2005, Miyaji et al. 2007). To this end we are developing the tools that will enable us 
to analyse large XMM data-sets, covering up to ~ 300 sq. degrees of sky. 

Although each of the previously discussed components of our project (Hubble relation 
using HH-like starburst galaxies and angular /spatial clustering of X-ray AGNs) will provide 
interesting and relatively stringent constraints on the cosmological parameters, especially under 
our anticipation that we will reduce significantly the corresponding random and systematic 
errors, it is the combined likelihood of these two type of analyses which will enable us to break 
the known degeneracies between cosmological parameters and determine with great accuracy 
the dark energy equation of state (see preliminary analysis of Basilakos & Plionis 2005; 2006). 

Below we present the basic methodology of each of the two main components of our proposal, 
necessary in order to constrain the dark energy equation of state. 

2. Cosmological Parameters from the Hubble Relation 

It is well known that in the matter dominated epoch the Hubble relation depends on the 
cosmological parameters via the following equation: 



which is simply derived from Friedman's equation. We remind the reader that Q m and Q w are 
the present fractional contributions to the total cosmic mass-energy density of the matter and 
dark energy source terms, respectively. 

Supernovae SNIa are considered standard candles at peak luminosity and therefore they have 
been used not only to determine the Hubble constant (at relatively low redshifts) but also to 
trace the curvature of the Hubble relation at high redshifts (see Riess et al. 1998, 2004, 2007; 
Perlmutter et al. 1998, 1999; Tonry et al. 2003; Astier et al. 2006; Wood-Vasey et al. 2007; 
Davis et al. 2007; Kowalski et al 2008). Practically one relates the distance modulus of the 
SNIa to its luminosity distance, through which the cosmological parameters enter: 



The main result of numerous studies using this procedure is that distant SNIa's are dimmer on 
average by 0.2 mag than what expected in an Einstein-deSitter model, which translates in them 
being ~ 10% further away than expected. 

The amazing consequence of these results is that they imply that we live in an accelerating 
phase of the expansion of the Universe, an assertion that needs to be scrutinized on all possible 
levels, one of which is to verify the accelerated expansion of the Universe using an alternative 
to SNIa population of extragalactic standard candles. Furthermore, the cause and rate of the 
acceleration is of paramount importance, ie., the dark energy equation of state is the next 
fundamental item to search for and to these directions we hope to contribute with our current 
project. 

2.1. Theoretical Expectations: 

To appreciate the magnitude of the Hubble relation variations due to the different dark energy 
equations of state, we plot in Figure 1 the relative deviations of the distance modulus, Afi, of 
different dark-energy models from a nominal standard (w = — 1) A-cosmology (with £l m = 0.27 
and £l\ = 0.73), with the relative deviations defined as: 




(3) 



fj, = m — M = 5 log d^ + 25 where 




A/X = fl\ — /imodel • 



(5) 
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Figure 1. Left Panel: The expected distance modulus difference between the dark-energy 
models shown and the reference A-model (w = —1) with Q m = 0.27. Right Panel: The expected 
distance modulus differences once the £l m -w(z) degeneracy is broken (imposing the same f2 m 
value as the comparison model). 



The parameters of the different models used are shown in Figure 1. As far as the dark- energy 
equation of state parameter is concerned, we present the deviations from the standard model 
of two models with a constant w value and of two models with an evolving equation of state 
parameter, utilizing the form of eq[2j In the left panel of Figure 1 we present results for selected 
values of Q m , while in the right panel we use the same dark-energy equations of state parameters 
but for the same value of f2 m (= 0.27) (ie., we eliminate the degeneracy). 
Three important observations should be made from Figure 1: 

(i) The relative magnitude deviations between the different dark-energy models are quite small 
(typically ^0.1 mag), which puts severe pressure on the necessary photometric accuracy 
of the relevant observations. 

(ii) The largest relative deviations of the distance moduli occur at redshifts z^, 1.5, and thus 
at quite larger redshifts than those currently traced by SN la, and 

(iii) There are strong degeneracies between the different cosmological models at redshifts z<^ 1, 
but in some occasions even up to much higher redshifts (one such example is shown in 
Fi gure 1 between the models with (£l m , wq,wi) = (0.31,-1,0) and (0.29,-1,0.2). 

Luckily, such degeneracies can be broken, as discussed already in the introduction, by using 
other cosmological tests (eg. the clustering of extragalactic sources, the CMB shift parameter, 
BAO's, etc). Indeed, current evidence overwhelmingly show that the total matter content of the 
universe is within the range: 0.2<; £l m ^> 0.3, a fact that reduces significantly the degeneracies 
between the cosmological parameters. 

2.2. Larger numbers or higher redshifts ? 

In order to define an efficient strategy to put stringent constraints on the dark- energy equation 
of state, we have decided to re-analyse two recently compiled SNIa samples, the Davies et al. 




Figure 2. Left Panel: SNIa distance moduli as a function of redshift. Right Panel: Distance 
moduli difference between the A-model and the SNIa data. The blue line is the corresponding 
difference between the reference (w = —1) and the w = —0.85 dark-energy models. 



(2007) [hereafter D07\ compilation of 192 SNIa (based on data from Wood-Vasey et al. 2007, 
Riess et al. 2007 and Astier et al. 2007) and the UNION compilation of 307 SNIa (Kowalski et 
al. 2008). Note that the two samples are not independent since most of the D07\s included in 
the UNION sample. 

Firstly, we present in the left panel of Figure 2 the UNION SNIa distance moduli overploted 
(red-line) with the theoretical expectation of a flat cosmology with ($l m ,w) = (0.27, —1). In the 
right panel we plot the distance moduli difference between the SNIa data and the previously 
mentioned model. To appreciate the level of accuracy needed in order to put constraints on the 
equation of state parameter, we also plot the distance moduli difference between the reference 
(tt m ,w) = (0.27, -1) and the (tt m ,w) = (0.27, -0.85) models (thin blue line). 

We proceed to analyse the SNIa data by defining the usual likelihood estimatoid as: 

£ SNIa (c)ocexp[- X | NIa (c)/2] (6) 
where c is a vector containing the cosmological parameters that we want to fit for, and 

y h (z. t ,c)-/i obs ( Zi ) 



XsNIa( c ) = 



(7) 



i=l 

where /i th is given by eq.(3), Zi is the observed redshift and <7j the observed distance modulus 
uncertainty. Here we will constrain our analysis within the framework of a flat (f2 tot = 1) 
cosmology and therefore the corresponding vector c is: c = [Q m , Wo, w±\. We will use only SNIa 
with z > 0.02 in order to avoid possible problems related to redshift uncertainties due to the 
local bulk flow. 

• Larger Numbers? The first issue that we wish to address is how better have we done in 
imposing cosmological constraints by increasing the available SNIa sample from 181 to 292 



2 Likelihoods are normalized to their maximum values. 
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Figure 3. Solution space of the fits to the Cosmological parameters using either of the two 
SNIa data sets. The contours are plotted where — 21n£/£ max is equal to 2.30, 6.16 and 11.83, 
corresponding to la, 2a and 3a confidence level. 



(excluding the z < 0.02 SNIa), ie., increasing the sample by more than 60%. In Table 
1 we present various solutions using each of the two previously mentioned samples. Note 
that since only the relative distances of the SNIa are accurate and not their absolute local 
calibration, we always marginalize with respect to the internally derived Hubble constant 
(note that fitting procedures exist which do not need to a priori marginalize over the 
internally estimated Hubble constant; eg., Wei 2008). 

We present our results in Figure 3. The left panel shows the internally derived Hubble 
constant for each of the to samples (and for w = — 1). It is evident that the derived 
values, which we use to free our analysis of the Hq dependence of the distance modulus, are 
well constrained (ho = 0.6515 ± 0.007 and 0.7015 ± 0.005 for the D07 and UNION sample, 
respectively). The right panel shows the main results of interest. The size of the welhknown 
banana shape region of the (Q m ,w) solution space is almost identical for both samples of 
SNIa. 

A first conclusion is therefore that the increase by ~ 60% of the UNION sample has not 
provided more stringent constraints to the cosmological parameters. Rather there appears 
to be an unexpected lateral shift of the contours towards higher values of fi m and lower 
values of w, within, however, the 1 a contour of the solution space. In order to verify that 
the larger number of SNIa's in the UNION sample are not preferentially located at low-z's 
- in which case we should have not expected more stringent cosmological constraints using 
the latter SNIa sample - we plot in Figure 4 the normalized redshift frequency distributions 
of the two samples (left panel) and the relative increase of SNIa's as a function of redshift 
in the UNION sample with respect to the D07 sample in the right panel. It is evident that 
the larger number of UNION SNIa's are distributed in all redshifts, except for z^ 1, where 
there is no appreciable increase of SNIa numbers. 

We already have a strong hint, from the previously presented comparison between the 
D07 and UNION results, that increasing the number of Hubble relation tracers, covering 
the same redshift range and with the current level of uncertainties as the available SNIa 




Table 1. Fitting the SNIa data (with z > 0.02 in order to avoid local bulk flow effects) to flat 
cosmologies. Note that for the case where c = (Q, m ,w) (ie. , last row), the errors shown are estimated 
after marginalizing with respect to the other fitted parameter. 





D07 UNION 


free params 


w n m Xmin/ df 


W ^rn. Xmin/ df 




-1 0.280+^ 187.03/180 
— 1.025±g;g|g 0.292 ±0.018 187.02/179 


-1 0.2801K 301.93/291 
-1.212 ±0.050 0.355 ±0.015 301.11/290 



samples, appears to be a futile avenue in constraining further the cosmological parameters. 

• Lower uncertainties or higher-z 's: We now resort to a Monte-Carlo procedure which will 
help us investigate which of the following two directions, which bracket many different 
possibilities, would provide more stringent cosmological constraints: 

— Reduce significantly the distance modulus uncertainties of SNIa, tracing however the 
same redshift range as the currently available samples, or 

— use tracers of the Hubble relation located at redshifts where the models show their 
largest relative differences (z<; 2), with distance modulus uncertainties comparable to 
that of the highest redshift SNIa's ((cr M ) ~ 0.4) 

The Monte-Carlo procedure is based on using the observed high-z SNIa distance modulus 
uncertainty distribution (c^) and a model to assign random /i-deviations from a reference 
H(z) function, that reproduces exactly the banana-shaped contours of the (p, m ,w) solution 
space of Figure 3 (right panel). Indeed, after a trial and error procedure we have found 
that by assigning to each UNION SNIa (using only their redshift) a distance modulus 
deviation from the reference model (in this case we use: £l m = 0.35, w = —1.21), 
using a Gaussian with zero mean and variance given by observed (a^) 2 ), and using as the 
relevant individual distance modulus uncertainty, which enters as a weight in the x 2 of eq.8, 
the following: af = y/ (1.15fii) 2 + <p 2 , with (j) a random Poisson deviate within [—0.01,0.01], 
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Figure 5. Left Panel: Comparison between the UNION SNIa constraints (black contours) and 
those derived by a Monte-Carlo procedure designed to closely reproduce them (for clarity we 
show only contours corresponding to 1 and 3 a confidence levels). Right Panel: The UNION 
SNIa distance modulus deviations from the best fit model - (£l m ,w) = (0.35,-1.21) - and a 
random realization of the model deviations (see text). 



we reproduce exactly the banana-shaped solution range of the reference model (details will 
appear elsewhere). This can be seen clearly in the upper left panel of Figure 5, where we 
plot the original UNIO iV solution space (black contours) and the model solution space (red 
contours). In the right panel we show the distribution of the true SNIa deviations from the 
best fitted model (Table 1) as well as a random realization of the model deviations. 
Armed with the above procedure we can now address the questions posed previously. Firstly, 
we reduce to half the random deviations of the SNIa distance moduli from the reference 
model (with the corresponding reduction of the relevant uncertainty, cjj). The results of the 
likelihood analysis can be seen in the left panel of Figure 6. There is a reduction of the range 
of the solution space, but indeed quite a small one. Secondly, we have increase artificially 
the high-z tracers by 88 objects (by additionally using the z > 0.65 SNIa and adding a 
6z = 2 to their redshift). Note that the new tracers are distributed between 2.65^ z<^ 3.55, 
ie., in a range where the largest deviations between the different cosmological models occur 
(see Figure 1). The deviations from the reference model of these additional SNIa are based 
on their original [i uncertainty distribution (ie., we have assumed that the new high-z 
tracers will have similar uncertainties as their z^, 0.65 counterparts, which is {a^) ~ 0.38). 
We now find a significantly reduced solution space (right panel of Figure 6), which shows 
that indeed by increasing the H(z) tracers by a few tens, at those redshifts where the 
largest deviations between models occur, can have a significant impact on the recovered 
cosmological parameter solution space. 

The main conclusion of the previous analysis is that the best strategy to decrease the 
uncertainties of the cosmological parameters based on the Hubble relation is to use standard 
candles which trace also the redshift range 2<^ z<? 4. Below we present such a possibility by 
suggesting an alternative to the SNIa standard candles, namely HH-like starburst galaxies (eg. 
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Figure 6. Comparison of the present SNIa constraints (red contours) with Left Panel: those 
derived by reducing to half their uncertainties (black contours) and Right Panel: with those 
derived by adding a sample of 88 high-z tracers (2.5<^ z<^ 3.6) with a distance modulus mean 
uncertainty of a a ~ 0.38 (black contours). 



Melnick 2003; Siegel et al. 2005). 

2.3. Hubble Relation using HH-like starburst galaxies 

We now reach our suggestion to use an alternative and potentially very powerful technique to 
estimate cosmological distances, which is the relation between the luminosity of the Hp line 
and the stellar velocity dispersion, measured from the line-widths, of HII regions and galaxies 
(Terlevich & Melnick 1981, Melnick, Terlevich & Moles 1988). The cosmological use of this 
distance indicator has been tested in Melnick, Terlevich & Terlevich (2000) and Siegel et al 
(2005) (see also the review by Melnick 2003). It is the presence of O and B-type stars in HII 
regions that causes the strong Balmer line emission, in both H a and Hp. Furthermore, the 
fact that the bolometric luminosities of HII galaxies are dominated by the starburst component 
implies that their luminosity per unit mass is very large, despite the fact that the galaxies are 
low-mass. Therefore they can be observed at very large redshifts, and this fact makes them 
cosmologically very interesting objects. Furthermore, it has been shown that the L(Hp) — a 
correlations holds at large redshifts (Koo et al. 1996, Pettini et al. 2001, Erb et al. 2003) and 
therefore they can be used to trace the Hubble relation at cosmologically interesting distances. 
One of the most important prerequisites in using such relations, as distance estimators, is the 
accurate determination of their zero-point. To this end, Melnick et al (1988) used giant HII 
regions in nearby late-type galaxies and derived the following empirical relation (using a Hubble 
constant of Hq = 71 km/sec/Mpc): 

log L(H p ) = log M z + 29.60 with M z = a 5 /{0/H) (8) 

where O/H is the metallicity. Based on the above relation and the work of Melnick, Terlevich 
& Terlevich (2000), the distance modulus of HII galaxies can be derived from: 



H = 2.51og(a 5 /F jff/3 ) - 2.51og(0/#) - A H p - 26.44 
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Figure 7. Comparison between the UNION SNIa constraints and those derived by using the 
15 high- z starburst galaxies of Siegel et al. (2005). Although, the latter constraints are weak 
and mostly inconsistent with the former, this plot serves to indicate the potential of using 
high-z starburst galaxies (once of course we have reduced significantly their distance modulus 
uncertainties). 

with Fjjp and Ajjp are the flux and extinction in Hp. The rms dispersion in distance modulus 
was found to be ~ 0.52 mag. The analysis of Melnick, Terlevich & Terlevich (2000) has 
shown that most of this dispersion (~ 0.3 mags) comes from observational errors in the 
stellar velocity dispersion measurements, from photometric errors and metallicity effects. It 
is therefore important to understand and correct the sources of random and systematic errors 
of the L(Hp) — a relation, and indeed with the availability of new observing techniques and 
instruments, we hope to reduce significantly the previously quoted rms scatter. 

A few words are also due to the possible systematic effects of the above relation. Such effects 
may be related to the age of the HII galaxy (this can be dealt with by putting a limit in the 
equivalent width of the Hp line, eg. EW(Hp) > 25 Angs; see Melnick 2003), to extinction, to 
different metallicities and environments. Also the EW(Hp) of HH-like galaxies at intermediate 
and high redshifts are smaller than the local ones, a fact which should be taken into account. 

We have commenced an investigation of all these effects by using high-resolution spectroscopy 
of a relatively large number of SDSS low-z HII galaxies, with a range of Hp equivalent widths, 
luminosities, metal content and local overdensity, in order to reduce the scatter of the HII- 
galaxy based distance estimator to about half its present value, ie., our target is ~ 0.25 mag. 
We will also define a medium and high redshift sample (see Pettini et al. 2001; Erb et al. 2003), 
consisting of a few hundred objects distributed in the different high-redshift bins, which will 
finally be used to define the high-z Hubble function. 

Summarizing, the use of HII galaxies to trace the Hubble relation, as an alternative to the 
traditionally used SN la, is based on the following facts: 

(a) local and high-z HH-like galaxies and HII regions are physically very similar systems 
(Melnick et al 1987) providing a phenomenological relation between the luminosity of the 
Hp line, the velocity dispersion and their metallicity as traced by O/H (Melnick, Terlevich 
& Moles 1988). Therefore HH-like starburst galaxies can be used as alternative standard 



candles (Melnick, Terlevich & Terlevich 2000, Melnick 2003; Siegel et al. 2005) 

(b) such galaxies can be readily observed at much larger redshifts than those currently probed by 
SNIa and 

(c) it is at such higher redshifts that the differences between the predictions of the different 
cosmological models appear more vividly. 

Already a sample of 15 such high-z starburst galaxies have been used by Siegel et al. (2005) 
in an attempt to constrain cosmological parameters but the constraints, although in the correct 
direction, are very weak. Here we perform our own re-analysis of this data-set and the resulting 
constraints on the £l m , w plane (for a flat geometry) can be seen in Figure 7. Note that imposing 
w = —1, our analysis of the Siegel et al (2005) data set provides O m = 0.10 ± 0.05, which is 
towards the lower side of the generally accepted values. Comparing these HH-based results to 
the present constraints of the latest SNIa data (D07 and UNION) clearly indicates the necessity 
to: 

• re-estimate carefully the local zero-point of the L(Hp) — a relation, 

• suppress the HH-galaxy distance modulus uncertainties, 

• increase the high-z starburst sample by a large fraction, 

• make sure to select high-z bona-fide HH-galaxies, excluding those that show indications of 
rotation (Melnick 2003). 

3. The Clustering of high-z X-ray AGN 

X-ray selected AGNs provide a relatively unbiased census of the AGN phenomenon, since 
obscured AGNs, largely missed in optical surveys, are included in such surveys. Furthermore, 
they can be detected out to high redshifts and thus trace the distant density fluctuations 
providing important constraints on super-massive black hole formation, the relation between 
AGN activity and Dark Matter (DM) halo hosts, the cosmic evolution of the AGN phenomenon 
(eg. Mo &; White 1996, Sheth et al. 2001), and on cosmological parameters and the dark-energy 
equation of state (eg. Basilakos & Plionis 2005; 2006). 

The earlier ROSAT-based analyses (eg. Boyle <fe Mo 1993; Vikhlinin & Forman 1995; Carrera 
et al. 1998; Akylas, Georgantopoulos, Plionis, 2000; Mullis et al. 2004) provided conflicting 
results on the nature and amplitude of high-z AGN clustering. With the advent of the XMM 
and Chandra X-ray observatories, many groups have attempted to settle this issue, but in vain. 
Different surveys have provided again a multitude of conflicting results, intensifying the debate 
(eg. Yang et al. 2003; Manners et al. 2003; Basilakos et al. 2004; Gilli et al. 2005; Basilakos et 
al 2005; Yang et al. 2006; Puccetti et al. 2006; Miyaji et al. 2007; Gandhi et al. 2006; Carrera 
et al. 2007). However, the recent indications of a flux-limit dependent clustering appears to 
remove most of the above inconsistencies (Plionis et al. 2008; see also Figure 8). 

Furthermore, there are indications for a quite large high-z AGN clustering length, reaching 
values ~ 15 — 18 h^ 1 Mpc at the brightest flux-limits (eg., Basilakos et al 2004; 2005, Puccetti 
et al. 2006, Plionis et al. 2008), a fact which, if verified, has important consequences for the 
AGN bias evolution and therefore for the evolution of the AGN phenomenon (eg. Miyaji et al. 
2007; Basilakos, Plionis &: Ragone-Figueroa 2008). An independent test of these results will 
be to establish that the environment of high-z AGN is associated with large DM haloes, which 
being massive should be more clustered. 

Below we justify the necessity for a large-area XMM survey in order to unambiguously 
determine the clustering pattern of high-z (z ~ 1) X-ray AGNs. We further show that such 
measurements can be used to put strong cosmological constraints (see for example Basilakos <fe 
Plionis 2005; 2006), and help break the tt m , w degeneracies. 
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Figure 8. The angular correlation scale, 0q, as a function of different survey characteristic flux, 
defined as that corresponding to half the respective survey area-curves. Most results appear 
to be consistent with the clustering flux-limit dependence, found from the CDF-N and CDF-S 
(from Plionis et al. 2008). 



3.1. Biases and Systematics: 

It is important to understand and overcome the shortcomings and problems that one is facing 
in order to reliably and unambiguously determine the clustering properties of the X-ray selected 
AGNs. Below we list the source of such problems, some of which can be rectified by considerably 
increasing the X-ray survey area. 

• Cosmic Variance: Is the volume surveyed large enough to smooth out inhomogeneities of the 
large-scale distribution of AGNs? (for example see Stewart et al. 2007). Closely related to 
this problem is the so-called integral constraint, which practically depends on the unknown 
true mean density of the cosmic sources under study. If the area is small enough, then the 
mean density, estimated from the survey itself, is way-out of its true value and thus the 
usual correlation function analysis will impose the observed mean number density as the 
true one (an example of this is the CDF-S were a large number of superclusters at z ~ 0.7 
are found; see Gilli et al. 2003). This usually results into an underestimation of the true 
correlation amplitude and a shallower zero-crossing of the estimated £(r) or w{0). A source 
of the observed scatter between the presently available surveys (see Figure 8) could well 
be the cosmic variance. These problems, however, are rectified with the large-area XMM 
survey proposed. 

• The amplification bias which can enhance artificially the clustering signal due to the 
detector's PSF smoothing of source pairs with intrinsically small angular separations (see 
Vikhlinin & Forman 1995; Basilakos et al. 2005). This problem can affect clustering 
results if at the median redshift of the sources under study the XMM PSF angular 
size corresponds to a rest-frame spatial scale comparable to the typical source pair-wise 
separations. Furthermore, the possible variability of the PSF size through-out the XMM 
fields can have an additional effect. This should be modeled and tested with Monte-Carlo 
simulations in order to establish the extent to which the clustering results are affected. In 



large-area surveys it is necessary to take good-care of this effect when using source pairs 
that belong to different XMM pointings. 

• Reliable production of random source catalogues: This is an issue which is extremely 
important and not appreciated at the necessary extent. The random catalogues, with 
which the observed source-pairs are compared, should be produced to account for all 
systematic effects from which the observations suffer, among which the different positional 
sensitivity and edge effects of each XMM pointing. Furthermore, a reliable log N — log S 
distribution (theoretically motivated or observationally determined) should be reproduced 
in the random "XMM pointings" and the random sources should be observed following the 
same procedure as in the true observations. Random positioned sources with fluxes lower 
than that corresponding to the particular position of the sensitivity map of the particular 
XMM pointing should be removed from each random catalogue. 

An optimal approach to unambiguously determine the clustering pattern of X-ray selected 
AGNs would be to determine both the angular and spatial clustering pattern. The reason being 
that various systematic effects or uncertainties enter differently in the two types of analyses. 
On the one side, using w{9) and its Limber inversion, one by-passes the effects of redshift-space 
distortions and uncertainties related to possible misidentification of the optical counter-parts 
of X-ray sources. On the other side, using spectroscopic or accurate photometric redshifts to 
measure £(r) or w p {6) one by-passes the inherent necessity, in Limber's inversion of w(9), of 
the source redshift-selection function (for the determination of which one uses the integrated 
X-ray source luminosity function, different models of which exist). For the inversion to work it 
is also necessary to model the spatial correlation function as a power law, to assume a clustering 
evolution model, which is taken usually to be that of constant clustering in comoving coordinates 
(eg. de Zotti et al. 1990; Kundic 1997) and also to assume a cosmological model. Limber's 
inversion then reads: 



where e = 7 — 3 for the constant clustering in comoving coordinates model, x(z) is the proper 
distance, E(z) = ^/f2 m (l + z) 3 + and = T(^)T( :1 ^-)/r(^). As noted previously, for the 
inversion to be possible it is necessary to know the X-ray source redshift distribution, dN/dz, 
and the total number, N, of the X-ray sources, which can be determined by integrating the X-ray 
source luminosity function above the minimum luminosity that corresponds to the particular 
flux-limit used. 

3.2. X-ray surveys 

In order to reach a flux-limit for which the soft-band clustering appears to converge to its final 
value (due to the flux-limit-clustering correlation, revealed in Plionis et al. 2008; see also Figure 
8) we plan to analyse all >10 ksec XMM pointings, which will allow us to reach a flux- limit of 
~ 2 x 10~ 15 erg/sec/cm 2 in the soft (0.5-2 keV) and ~ 10~ 14 erg/sec/cm 2 in the hard (2-10 keV) 
bands, respectively. Such an exposure time will finally provide (taking in to account realistci 
observational effects) ~250 soft-band and ~100 hard-band X-ray sources per deg 2 , according 
to the Kim et al. (2007) log N — log S and since all publicly available XMM pointings, with 
exposure time more than 10 ksec, add to ~ 300 non-contiguous sq. degrees, implies a resulting 
sample of ~ 75000 soft and ~ 30000 hard X-ray sources. These numbers will allow us to derive 
with great accuracy the small-separation angular correlation function, the Limber's inversion of 
which can provide their spatial correlation function. Furthermore, around 100 non-contiguous 
sq.degrees of the previous survey (except for a few contiguous regions with areas between 2 and 
10 sq.degrees each) are covered also by the SDSS, providing crude photo-zs. This will allow 




(10) 



us to derive the angular correlation function in distinct redshift bins in the range 0.5 < z < 2 
and thus quantify the evolution of the bias of X-ray selected AGN, an important ingredient in 
disentangling the cosmological parameters. Note finally that ~50 sq. degrees are covered also 
by the publicly available UKIDSS Qhttp : //www . ukidss . o rg/) which provide 3 near-IR colours 
and thus for a subsample of the previous data we will have relatively more accurate photo-zs 
allowing us to attempt to derive directly the spatial correlation function. 

Summarizing, the analysis of such large XXM surveys will allow us to unambiguously 
determine the soft and hard-band X-ray AGN clustering pattern, minimizing the biases and 
systematic effects discussed in the previous section, as well as to study the evolution of the AGN 
correlation function (utlizing photo-z's and dividing the angular sample into 2-3 z-bins). 

Finally, together with a large European consortium, we are planning to survey a 2-3 
contiguous sky areas adding to ~ 50 deg 2 and covered also by a large number of other 
multiwavelength surveys (like UKIDSS, NEWFIRM, etc), which will probably allow us to 
sample the long wavelength regime of the AGN correlation function, and thus measure Baryonic 
Accoustic Oscillations, which are extremely important for Dark-energy investigations (Peacock 
et al 2006; Albrecht et al. 2006). To obtain the necessary high photo- z accuracy, one may 
envision a wide-field multi-filter (say using ~20 broad and narrow-band filters) survey, specially 
tuned in order to obtain relatively high-accuracy photo-z's of z ~ 1 AGN. Thoughts for the 
aquisition of such an instrument, to be mounted on the 2.3m Hellenic Aristarchos telescope, are 
already been discussed in the Institute of Astronomy & Astrophysics of the National Observatory 
of Athens. 



3.3. Cosmological Parameter constraints from X-ray AGN Clustering 

The unambiguous determination of the correlation function of the z ~ 1 X-ray AGNs, even in 
angular space, will allow us to estimate with good precision (a) their relation to the underlying 
matter fluctuations at z ~ 1 (ie., their bias), (b) the evolution of their bias and therefore the 
mass of the DM haloes which they inhabit (eg., Miyaji et al. 2007; Basilakos, Plionis & Ragone- 
Figueroa 2008) and (c) put strong cosmological constraints on the f2 m , h or Q m , planes, while 
with the help of the high-z Hll-based Hubble relation on the Q m ,w(z) space. 

It is well known (Kaiser 1984; Benson et al. 2000) that according to linear biasing the 
correlation function of the AGN (or any mass-tracer) (£agn) and dark-matter one (£dm)> ar e 
related by: 

£AGx(r,z) = b 2 (z)£, DM (r,z) , (11) 

where b{z) is the bias evolution function (eg. Mo & White 1996, Matarrese et al. 1997, Basilakos 
& Plionis 2001; 2003; Basilakos, Plionis & Ragone-Figueroa 2008) 

A first outcome of the proposed correlation function analysis will be the accurate 
determination of the AGN bias at their median redshift (in our case z — 1), utilizing: 

/ \ 7/2 

b(z) = — D 3+e (z) with 7 = 1.8 and e = -1.2 , 
\ r o,m / 

where ro and ro, m are the measured AGN and dark-matter (from the P(k)) clustering lengths, 
respectively, while D{z) is the perturbation's linear growing mode. Then using a bias evolution 
model, one will be able to determine the mass of the DM halo within which such AGN live (eg. 
see Figure 9). 

Furthermore, we will compare the observed AGN clustering with the predicted, for different 
cosmological models, correlation function of the underlying mass, ^M( r , z )- To this end we can 
use the Fourier transform of the spatial power spectrum P(k): 

&mM) .(I±|^^p W ^, (12) 
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Figure 9. Comparison of the Basilakos et al (2008) b(z) evolution model with different 
observational data. Different line types represent different halo masses. Left Panel: optical 
galaxies (open points) with solid line corresponding to Mdm ~ 10 12 h^ 1 Mq, Lyman break 
galaxies (solid triangle), EROs (star), DRGs (open triangle) and 2dF radio sources (filled 
pentagon). The dot-dashed line corresponds to 7.7 xlO 13 h~ l Mq. Right Panel: optically 
selected quasars (open points and crosses), soft and hard X-ray point sources (open squares and 
solid diamonds; the large scatter corresponds to the uncertainty of the present day clustering 
results; see Plionis et al. 2008 and references therein). In the insert we plot, as solid points, the 
high-z SSRS DR5 QSOs and the same b(z) model that fits their lower redshift counterparts (ie., 
M DM = 10 13 h~ l Mq). 




Figure 10. Likelihood contours from the X-ray AGN clustering analysis of Basilakos &: Plionis 
(2005) in the following planes: (r2 m ,o"s) (right panel), (as, bo) (left panel). The contours plotted 
correspond to la, 2a and 3a confidence level. 




Figure 11. Left Panel: Likelihood contours on the £l m ,w plane from the X-ray AGN clustering 
analysis of Basilakos & Plionis (2005; 2006) (red contours) and the UNION SNIa analysis (black 
contours). Right Panel: The Joint likelihood contours. 

where k is the comoving wavenumber, P(k) = P k n T 2 (k) the CDM power-spectrum with scale- 
invariant (n = 1) primeval inflationary fluctuations, while the transfer function parameterization 
is as in Bardeen et al. (1986), with the corrections given approximately by Sugiyama (1995). 

Basilakos & Plionis (2005; 2006) have already used succesfully a standard maximum likelihood 
procedure to compare the measured XMM source angular correlation function from a relatively 
small (~ 2 sq.degrees survey; Basilakos et al. 2005) with the prediction of different spatially flat 
cosmological models, and derived interesting cosmological contraints (for flat and constant-u; 
cosmologies) . 

In Figure 10 we present the constraints provided from the preliminary analysis by Basilakos 
& Plionis (2006) of a ~ 2 sq.degrees XMM survey, on the present bias- factor (&o) of the X-ray 
AGN, on the normalization of the power spectrum (as) and on f2 m (marginalizing over different 
parameters). In Figure 11 (left panel) we present the (O m ,w) constraints provided by the X-ray 
AGN clustering analysis (red contours), once we have marginalized over the ag(~ 0.8) and the 
bias factor at the present time (~ 2). 

4. Joint Hubble-relation and Clustering analysis 

It is evident from Figure 11 (left panel) that w is degenerate with respect to Q m and that all 
the values in the interval —2 < w < —0.35 are acceptable within the la uncertainty. However, 
we break this degeneracy by adding the constraints provided by the Hubble relation technique, 
using here the UNION SNIa sample. We therefore perform a joint likelihood analysis, assuming 
that the two data sets are independent (which indeed they are) and thus we can write the joint 
likelihood as the product of the two individual ones. 

Our current joint likelihood analysis, once we impose h = 0.72 and a% = 0.8, provides quite 



stringent constraints of the Vt m and w parameters: 

n m = 0.28jl8;82 and id = -1.0 ± 0.1 . 

However, the uncertainty on w is still quite large, while the necessity to impose constraints on 
a more general, time-evolving, dark-energy equation of state (eq. 2) implies that there is ample 
space for great improvment and indeed the aim of the project, detailed in these proceedings, is 
exactly in this direction. 
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